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Abstract. 

We report the first measurements and detailed analysis of extreme ultraviolet 
(EUV) spectra (4 nm to 20 nm) of highly-charged tungsten ions W^'*+ to W^^+ 
obtained with an electron beam ion trap (EBIT). CoUisional-radiative modelling is 
used to identify strong electric-dipole and magnetic-dipole transitions in all ionization 
stages. These lines can be used for impurity transport studies and temperature 
diagnostics in fusion reactors, such as ITER. Identifications of prominent lines from 
several W ions were confirmed by measurement of isoelectronic EUV spectra of Hf, Ta, 
and Au. We also discuss the importance of charge exchange recombination for correct 
description of ionization balance in the EBIT plasma. 
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Future fusion reactors, such as ITER, are to reach temperatures of about 20 keV 
to 25 keV [1]. While hght elements (e.g., D, T, and C) will be completely ionized in 
the plasma core interior, heavy impurity ions will still possess a number of electrons. 
This could result in strong line emission, primarily in the x-ray region. Such radiation 
is the source of power losses that represent one of the major concerns for sustainable 
fusion. On the other hand, high resolution x-ray spectroscopy of impurities offers reliable 
measurements of important plasma parameters such as ion temperature Tj, rotation 
velocities in toroidal and poloidal directions, and electron temperature Tg [2]. Among 
all possible impurities in ITER, tungsten is expected to be the most abundant, since 
according to current planning the front surface of the divertor will be made of this 
element. 

Extreme ultraviolet (EUV) line emission of impurities has attracted less attention 
than the x-ray lines, in part because the flux of emitted photons in the EUV region is 
typically smaller due to lower transition probabilities. Feldman et al [3] recently pointed 
out that a number of EUV lines from highly-charged tungsten (ion charge z > 53) may 
be reliably recorded in the ITER plasma and used for diagnostics of plasma temperature 
and impurities transport. They proposed, in addition to the traditional grazing incidence 
spectrometer, a system of segmented multilayer-coated telescopes for registration of 
EUV lines, similar to what is used in solar corona diagnostics (see, e.g., [1]). Other 
possibilities for registration of EUV lines in ITER are also discussed in the literature 

The objective of the present work is to investigate the EUV spectra from those 
highly-charged W ions that will be abundant in the plasma of ITER. The measurements 
of the spectra were performed with the Electron Beam Ion Trap (EBIT) at the National 
Institute of Standards and Technology (NIST). EBIT is a versatile light source, capable 
of producing nearly any ion charge state of nearly any element. A fine control of the 
charge state distribution in the trap is due to a very narrow electron energy distribution 
function (EEDF) of the beam (width < 60 eV) j5j- A detailed description of the 
NIST EBIT can be found elsewhere [6]. The ions under study were injected into the 
EBIT using a multi-cathode, metal vapor vacuum arc (MEVVA) |7j designed for a 
rapid change of the injected element without downtime or alteration of experimental 
conditions. The EUV spectra between 4 nm and 20 nm were recorded with a grazing- 
incidence spectrometer described in detail in [8]. The instrument's resolution is about 
350, corresponding to a resolving limit of about 0.03 nm. The W spectra were calibrated 
with lines of highly- charged ions of Fe, with wavelengths taken from the NIST Atomic 
Spectra Database [9]. 

The measured spectra for six nominal electron beam energies Eb varying from 
8.8 keV to 25 keV are presented in figure [H A typical value of the beam current was 
about 150 mA. Unlike our recent work on low-energy EUV tungsten spectra [in], here 
we did not use a zirconium foil to filter out light from wavelengths above 25 nm; this 
resulted in significantly higher signal-to-noise ratio near the edges of the spectral range. 
For easier visual identification of the second-order lines, the shifted spectra (red line) 
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show the same experimental spectra with wavelengths multiplied by 2. The spectra also 
show a number of impurity lines, mainly from oxygen. The oxygen lines at 13.3 nm, 
15.0 nm, 17.2 nm, and 17.3 nm are marked by crosses in the spectrum for Eb = 9.3 keV 
(figure [1]). 

As with our previous studies of x-ray and EUV spectra of tungsten ions with 
lower charges {z = 37 to 47) [lOl [11], identification of the spectral lines relies upon 
coUisional-radiative modelling performed with the non-Maxwellian code NOMAD |12j . 
The Flexible Atomic Code (FAC) [I3], based on the relativistic model-potential method 
and jj-coupling scheme, was used to generate all relevant atomic data such as energy 
levels, wavelengths, transitions probabilities, and coUisional cross sections for the ions 
W^^^ through W^^+. Previous comparisons with the KBIT tungsten spectra at lower 
beam energies [lOl [TT] showed that FAC provides sufficiently accurate results for highly- 
charged ions. Using the FAC data, NOMAD calculates the ionization balance, level 
populations, and line intensities for the EBIT conditions. The EEDF in our simulations 
was modelled as a Gaussian peak with a full width at half maximum of 60 eV, and the 
electron density was set at = 2x10^^ cm~'^. The calculated spectra were convolved 
with a Gaussian instrumental function and relative efficiency curve of the spectrometer 

The identifications of the spectral lines, the measured and calculated wavelengths, 
transition types, and calculated radiative decay rates are presented in table [U The total 
uncertainties in the wavelengths reported in the table are typically dominated by the 
uncertainty in the spectrometer calibration, rounded up to one significant digit. An 
additional statistical uncertainty due to fitting the reported line centers has been added 
in quadrature. To facilitate classification, for each state we also report a calculated 
level number in the energy-ordered list of levels within the corresponding ion stage. All 
identified spectral lines correspond to An = transitions within the n = 3 shell. The 
LS'-coupling identifications of the upper and lower levels that were obtained using the 
Cowan code [l4j are only approximate due to strong spin-orbit interaction. 

A straightforward method to match experimental and calculated line intensities 
consists in varying the theoretical beam energy until a reasonable agreement is achieved. 
While the calculated wavelengths agree well with the measurements, our modelling 
required significantly lower beam energies to reproduce the observed relative line 
intensities. For instance, for the highest nominal beam energy of E^ = 25 keV, the 
best-fit theoretical energy was only 10.5 keV, while for E^ = 8.8 keV the fit energy 
was about 6.2 keV. Although space-charge effects in the trap may effectively reduce the 
beam energy, for the present values of the nominal voltage and beam current this would 
result in a less than 300 eV shift. 

Calculation of ionization balance in an EBIT plasma was a subject of several studies 
(e.g., fi5[ [T6l IT71 IT8] and references therein). Each ion is normally modelled as one 
state with no internal structure, and approximate formulas, such as Lotz formula for 
ionization cross sections, are used for description of atomic processes. In constrast, our 
coUisional-radiative modelling emphasizes the atomic physics component of simulations 
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Figure 1. (colour online) Experimental spectra at six beam energies from 8.8 keV to 
25 keV. For better identification, the thin red line shows the second order spectrum. 
Oxygen impurity lines are marked by crosses. 
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Table 1. Identified lines of highly- charged ions of tungsten. The experimental 
uncertainties for wavelengths are given in parentheses. The level number within the 
corresponding ion stage is shown in square brackets. 
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by implementing detailed descriptions of charge stages (about 1000 levels per ion) and 
relevant physical processes affecting populations of atomic levels. This is a necessary 
condition for the accurate modelling of well-resolved EBIT spectra. 

It is known (e.g., [10]) that the relative line intensities within a particular ionization 
stage are only weakly sensitive to variations in the EBIT beam energy. Therefore the 
observed disagreement with the simulations is most likely due to modification of the 
ionization distribution in the EBIT plasma. This discrepancy can be explained by charge 
exchange (CX) with neutral particles (mostly nitrogen and oxygen) that are always 
present in the trap. First, note that the only important electron-ion recombination 
process in the EBIT is the radiative recombination. Dielectronic recombination (DR) 
would generally be very important for the high ion charges in question; however, it is 
an energy-selective process, and for all cases considered here the electron beam energy 
is outside the range of importance for DR. Our present calculations show that the 
typical radiative recombination cross sections for the relevant ions of tungsten are on 
the order of 10~^^ cm^ (see also [19]). It is well known that the CX cross sections are 
normally much larger. The Classical Trajectory Monte-Carlo (CTMC) theory [20] gives 
the approximate dependence as: 

acx ^ z X 10"^^ cm^, (1) 

and thus for 60 to 70 times ionized atoms the CX cross section is close to 10~^^ cm^. 
This simple formula agrees within 40% with the semiempirical scaling law for charge 
exchange [TH [21] that is often used in charge distribution simulations of EBIT plasmas. 

The CX rate depends on the density of the neutrals A^^, relative velocity of neutrals 
and tungsten ions Vr, and the cross section acx- For simplicity, the product {(icxVr), 
averaged over the relative energy distribution of neutrals and W ions, can be replaced 
by the product of acx and effective averaged velocity Vr, so that the CX rate becomes 
Rex = N^acx^r- Both Nn and Vr are unknown parameters that may actually vary 
from measurement to measurement. (A slight enhancement of lower ionization stages at 
Eb = 9.0 keV vs. Eb = 8.8 keV may be due to this uncertainty.) However, since they 
enter only as a product, one can use this single parameter to fit the ionization balance. 
As there exist no calculations of CX for such highly-charged ions, we used the CX cross 
section of Eq. and varied the NnVr product to reach agreement with experimental 
spectra. The example of correspondence for the most complex spectrum at Eb = 8.8 
keV is presented in Fig. [2l One can see that the calculated spectrum that uses only one 
fit parameter (the theoretical beam energy was fixed at 8.8 keV) correctly reproduces 
relative intensities for practically all observed lines. The derived value of NnVr ~ 4x10^^ 
cm~^s~^ is consistent with general estimates for and Vr expected inside the EBIT. 
The best theoretical fit for the experimental spectrum at Eb = 25 keV is obtained with 
NnVr ~ 1x10-^^ cm'^s"-*^. 

As mentioned above, our simulations were performed at Ne = 2x10^^ cm^^. 
Without CX, the ionization balance is practically independent of Ne for a several-orders- 
of-magnitude span of densities. With CX, variations in Ng would certainly modify 
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Figure 2. (colour online) Comparison of experimental (top) and theoretical (bottom) 
spectra at the beam energy Eb — 8.8 keV. The second-order spectrum is shown by 
the red line. The positions of the strongest impurity lines from highly-ionized oxygen 
are indicated by dashed lines. 



ionization rates. However, since NnVr is a free parameter in our simulations, it can 
be appropriately adjusted to restore agreement with the experimental spectra. Thus, 
dependence on Nf, is effectively removed from the modelling. 

At the lowest energies Eb = 8.8 keV and 9.0 keV (see figure [1]), one can clearly 
identify two groups of strong lines, namely, in the range of 4 nm to 8 nm and 14.5 nm 
to 17.5 nm. These lines originate from ions W^^^ through W^°"*". As the beam energy 
increases, the ionization distribution shifts to higher charge states, and the three lines 
from Al-like W^^+ (A = 7.404 nm), Mg-like W^^+ (A = 7.991 nm), and Na-like W^3+ 
(A = 7.769 nm) ions become the most prominent. The variation of their intensities 
agrees with the changes in the beam energy. These three lines will be easily resolved 
in plasmas of fusion reactors, and in addition, they cover a range of only 0.6 nm, so 
that determination of their relative intensities will be less dependent on the spectral 
sensitivity of EUV spectrometers. 

To confirm our identification of these strong lines in W, we also measured the EUV 
spectra for Hf (Z=72), Ta {Z=73), and Au {Z=79). The nominal beam energy for Au 
was 14 keV, while the Hf and Ta spectra were recorded at Eb = 12 keV. This selection 
of energies provided ionization distributions similar to the 12 keV measurements for W. 
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Table 2. Measured and calculated wavelengths (in nm) of the strongest lines from A1-, 
Mg-, and Na-likc ions of Hf, Ta, and Au. The experimental wavelength uncertainty is 
± 0.002 nm for all lines in this table. 



Sequence 


Transition 




Hf 




Ta 


Au 


expt 


calc 


expt 


calc 


expt calc 


Al 


Ss^Sp - 3s3p2 4p,/2 


7.741 


7.7297 


7.570 


7.5598 


6.644 6.6459 


Mg 


3s2 ^So - 3s3p 


8.352 


8.3456 


8.169 


8.1652 


7.181 7.1902 


Na 


3s 2^1/2 - 3p2Pi/2 


8.107 


8.0973 


7.937 


7.9270 


6.999 7.0030 



Figure [3] shows the recorded spectra between 6 nm and 9 nm, including the W results. 
The measured and calculated wavelengths for Hf , Ta, and Au are presented in Table |2j 
As these lines correspond to the An = transitions, the leading term in the transition 
energy difference should be proportional to the net charge of the atomic core Zc = 
z + 1. This dependence is well confirmed in Fig. H] which shows the linear fit results 
for wavenumbers a. The fit a [cm~^] = A ■ - B agrees with the experimental data 
to within 0.1 % for the following fit coefficients (A,B): (3.056 xlO^ 5.428x10^) for the 
Al-like line, (2.791xl0^ 5.052x10^) for the Mg-like line, and (2.801xl0^ 5.043x10^) 
for the Na-like line. The derived coefficients can be used to predict wavelengths of these 
strong lines in other elements. 

There exist only a few theoretical works on calculated energy levels, wavelengths, 
and transition probabilities for W ions with 54 < z < 63, with most efforts related to 
Na-like W^^"*". Ali and Kim [23] used multiconfiguration Dirac-Fock (MCDF) method to 
calculate the wavelength of the 3s^3d "^0^/2 ^ 3s^3d ^.05/2 transition in W^^^, and their 
value of A = 15.959 nm is within the experimental uncertainty. The wavelength of the 
3s2 ^5*0 - 3s3p transition in Mg-like W^2+ 

was calculated at A = 8.0509 nm using the 
relativistic model-potential method [25j, compared to our measured value of 7.991(2) 
nm. Among the recent calculations for the 3si/2 - 3pi/2 transition in the Na-like ion, 
the MCDF [251 126] and many-body perturbation theory (MBPT) [27] results should be 
mentioned. The latter give a wavelength A = 7.7659 nm in very good agreement with the 
present experiment, and the MBPT transition probability A = 4.296 x 10^° s~^ agrees 
with our calculations within 2 %. The ab initio MCDF wavelength of Seely et al [25] A 
= 7.7479 nm is also close to the measured value, while Baik et al 's result of A = 7.8184 
nm [26] is higher by 0.05 nm. Finally, the wavelength predicted in [25] by extrapolation 
along isoelectronic sequence, A = 7.7661 nm, agrees well with the experiment. 

The low electron density in the EBIT, two to three orders of magnitude lower 
than that in typical tokamaks, allows one to observe weak spectral lines that might be 
collisionally quenched in fusion reactors. A recently discussed example is provided by the 
magnetic-octupole line in the Ni-like W^^^ [22], with a radiative transition probability 
on the order of 10^ s~^, that becomes quenched at densities above 10^^ cm~^. The 
higher degree of ionization for the lines considered here should be more favorable for their 
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Figure 3. (colour online) Experimental spectra for Hf (beam energy 12 keV), Ta (12 
keV), W (12 keV), and Au (14 keV). The three strongest lines originate from (short 
to long) A1-, Na-, and Mg-like ions, respectively. 



survival as the most important An=0 collisional cross sections decrease with an increase 
of ion charge. To test this assumption, we used the NOMAD code to solve the steady- 
state system of rate equations for approximately 6000 levels of the ions W^^"*" through 
W^^^ at an electron temperature Tg = 10^ eV and density = 10^^ cm~^. It was found 
that for all lines in table [1], including the weakest lines with transition probabilities A < 
10® s~\ the radiative decays are still the dominant depopulation process for the upper 
levels of corresponding transitions. This result remains valid for a large range of electron 
temperatures since the typical energy difference for collisional processes is much smaller 
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Figure 4. (colour online) Zc-dependence of wavenumbers (in cm"^) for the three 
strongest lines of Fig. [31 Al-like - circles, Mg-like - squares, Na-like - diamonds. The 
linear fits are shown by dashed lines. The error bars for wavenumbers are completely 
within the symbols. 



than Tg. Therefore, these hnes will not be coUisionally quenched in the ITER plasma. 

To conclude, we presented here the first measurements of EUV spectra from highly- 
charged [z = 54 to 63) ions of W that are important for plasma diagnostics in ITER- 
class fusion reactors. Using advanced coUisional-radiative modelling of line spectra from 
the non-Maxwellian plasma of EBIT, we identified 26 new electric-dipole and magnetic- 
dipole lines from 10 consecutive ionization stages W^^^ through W^'^"''. We also measured 
wavelengths of three important spectral lines for A1-, Mg-, and Na-like ions of Hf, Ta, 
and Au and used them to confirm our identifications of corresponding lines in W. These 
prominent lines may be used as tracer identifications in the ITER plasma interior. The 
increasing importance of charge exchange for ionization balance of high-Z ions in the 
EBIT plasmas was pointed out as well. A detailed description of the modelling and a 
discussion of various spectral features will be provided elsewhere. 
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